Can the Mott Insulator TiOCl be Metallized by Doping? 
A First-Principles Study. 
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We investigate the effect of Na intercalation in the layered Mott insulator TiOCl within the 
framework of density functional theory. We show that the system remains always insulating for all 
studied Na concentrations, and the evolution of the spectral weight upon Na doping is consistent 
with recent photoemission experiments. We predict the Na-doped superlattice structures, and show 
that substitutions of O by F, CI by S, or Ti by V (or Sc), respectively, fail to metallize the system. We 
propose a description in terms of a multiorbital ionic Hubbard model in a quasi-two-dimensional 
lattice and discuss the nature of the insulating state under doping. Finally, a likely route for 
metallizing TiOCl by doping is proposed. 

PACS numbers: 71.15.Mb,71.15.Pd,71.27.+a,71.30.+h,71.70.Ch 
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A successful route to drive a layered compound into 
a (usually unconventional) superconducting state is by 
_doping the system vifith additional electrons or holes [l|- 
Jn contrast to the cuprates 0, Q , Fe- based systems 
and cobaltates where 3d ions with occupations 
larger than 5 electrons are responsible for superconduc- 
tivity, in the layered nitride superconductors Q, Zr and 
Hf are in 4d^ and 5d^ configurations. In this work we in- 
vestigate whether unconventional superconductivity can 
be realized in early 3d layered compounds upon dop- 
ing. A good candidate is TiOCl which is, similar to 
the cuprates, a layered Mott insulator at room temper- 
ature 1] but with Ti'^^ in its 3d^ configuration; Ti is 
in the same column of the periodic table as Zr and Hf. 
According to these features, a fundamental question is 
whether TiOCl can be a metal or even superconductor 
by doping. 

In fact, the undoped TiOCl, formed by stacked bilayers 
of TiO along the c axis separated by CI ions, has already 
attracted considerable interest |9l4l3| due to the appear- 
ance of a zero-field structural incommensurate state. Fur- 
thermore, great effort has been put into metalli zing the 
undoped TiOCl by applying hydrostatic pressure 14. isj. 

Very recently, samples of TiOCl with Na intercalated 
into the interstices of the TiO bilayers were prepared in 
order to drive the system into a metallic or even su- 
perconducting state upon electron doping |l6|. How- 
ever, photoemission experiments (PES) show that a pro- 
nounced gap remains at the Fermi level for all concentra- 
tions of Na, indicating that the insulating state always 
survives electron doping [l6| . In view of these results, it is 
essential to understand the microscopic origin of this puz- 
zling insulating state upon doping and, thereafter, find 
a more promising recipe for metallizing TiOCl. More- 
over, considerable interest in insulator-insulator transi- 
tions [l3| motivates us to investigate the nature of the 
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FIG. 1: (Color online) Superlattices at 6 different concentra- 
tions of Na doping, (a) undoped case. The atoms where sub- 
stitutions of CI by S, O by F and Ti by V(Sc) are performed 
are also specified, (b) 12.5 % Na doping: Naj/gTiOCl. The 
8 Ti ions in the cell are labelled by numbers 1 to 8. (c) 25 % 
doping: Nai/4TiOCl, (d) 50 % doping: Naj/zTiOCl, (e) 75 
% doping: Na3/4TiOCl, and (f ) 100 % doping: NaTiOCl. 



doped insulating state. 

In this Letter, we performed ab initio molecular dy- 
namics studies with the Car-Parrinello projector- 
augmented wave (CP-PAW) method [l3| to study the 
effect of Na do ping on TiOCl within the Parrinello- 
Rahman scheme [20| . The Perdew-Burke-Ernzerhof gen- 
eralized gradient approximation (GGA)-|-[/ to the den- 
sity functional theory (DFT) is employed in order to have 
the correct insula ting behavior and lattice structure in 
the undoped case [3 [ll, [ill-ii]. 

Our findings can be summarized as follows: For all 
Na concentrations studied, the unusual insulating state 
persists and the change of spectral weight at elevated 
Na concentration agrees well with the PES results flit . 
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FIG. 2: (Color online) Total density of states (DOS) at the 
6 different concentrations of Na doping shown in Fig. [T] (a) 
undoped case: TiOCl, (b) Nai/gTiOCl, (c) Nai/4TiOCl, (d) 
Nai/aTiOCl, (e) Na3/4TiOCl, and (f) NaTiOCl. 



Analysis of this insulating state shows that the apro- 
priate microscopic model is no longer a dominantly 
quasi-one-dimensional single-band model valid in the un- 
doped case [i, 11- 13 1 but a two-dimensional multi-orbital 
model. The effect of Na doping is twofold; (i) it strongly 
modifies the crystal field splitting of Ti d states by the 
Na+ ion and (ii) it induces a coexistence of Ti'^+ and Ti^+ 
ions which is further stabilized by long-range Coulomb 
interaction as usually observed in mixed-valence com- 
24| . Together with the on-site Coulomb repul- 



pounds 

sion as well as the Hund's rule coupling, these effects play 
a crucial role in forming the doped symmetry-breaking 
insulating phase. We further predict the Na-doped super- 
lattice structures, and we also find that doping TiOCl by 
substituting O by F, CI by S, or Ti by V (or Sc) does not 
result in a metal. While the insulating state of the F/S 
substitution case can be ascribed to the same microscopic 
origin as in the Na-doped case, the long-range Coulomb 
interactions play a minor role in the V/Sc substitution. 

For our calculations we used supercells obtained by 
doubling the primitive unit cell (space group Pmmn, 
two formula units per unit cell) in both a and b direc- 
tions. The molecular formula under doping is defined as 
Naa;Tii_yVy(Scy)0i_2F2Cli_u,Su,, where x, y, z, and w 
are the doping concentrations (Fig.[T](a)). We performed 
carefully converged relaxations of lattice parameters and 
atomic positions for each doping case with high energy 
cutoffs of 45 Ry and 180 Ry for the wave functions and 
charge density expansion, respectively. 

We first study the case of Na^: TiOCl for six different 
values a:=0, 1/8, 1/4, 1/2, 3/4, 1. Since no experimental 
data for the Na positions are available, numerous possi- 
bilities have to be tested by ab initio molecular dynamics. 
However, according to Pauling's rule that the coordina- 
tion number is determined by the radius ratio (25| . the 
most probable Na positions are in the cages of 5 CI and 
1 O as shown in Figs.[T] (b)-(f). In fact, by performing ab 
initio molecular dynamics starting from several different 
initial Na positions, we find that the Na atoms always 
fall into this cage to reach the coordination number of 6. 
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FIG. 3: (Color online) Partial DOS for (a), (d) 12.5% 
Na doping: Nai/gTiOCl, (b), (e) 12.5% F substitu- 
tion: TiOr/gFi/gCl, and (c), (f) 12.5% S substitution: 
TiOClr/sSi/g. In (a), (b), and (c) the total Ti-resolved DOS 
are shown. The grey (orange) curve denotes the DOS for the 
Ti ion nearest to the doping ion (Ti2 for the Na- and F-doped 
cases, Ti3 for the S-doped case). The black (blue) curve de- 
notes the DOS for the rest of Ti ions. In (d), (e), and (f) the 
corresponding orbital-resolved DOS are shown. Below Ef the 
contributions of Ti2 for the electron-doped cases are explicitly 
marked. 



The remaining uncertainty is given by the various com- 
binations of locating Na atoms ipto the 8 cages of the su- 
percell. We performed ab initio molecular dynamics for 
all the different combinations at various doping concen- 
trations and compared the total energies. Fig. [T] presents 
the final stable lattice structures with lowest total ener- 
gies for the six Na concentration values x studied [iij . 
The supercell structures as a function of Na doping are 
shown in the figure and call for further experimental X- 
ray diffraction studies. 

In Fig. [5] we present the total density of states (DOS) 
for the six Na-doped TiOCl cases considered. The DOS 
of the undoped case (Fig. [2] (a)) is correctly reproduced 
with an energy gap around 2 eV, and the peak close to 
the Fermi level is dominantly of dx^-y^ character. Here 
we chose the local coordinate frame as cc || 6, y \\ c and 
z II a with d.j.2^y2 , dxz, dyz forming the t2g bands and dxy, 
dz2 the eg bands. With Na doping (Figs. [2] (b)-(f)) the 
whole spectral weight is suddenly shifted towards lower 
energies, and an additional peak with weight 2x appears 
close to the Fermi level as compared to the undoped case. 
The separation of the two Ti peaks is about 1 eV for 
a; = 1/8 and a gap always persists at the Fermi level in 
all Na-doped cases. All the findings are consistent with 
the experimental observations [16i |. 

In order to understand the microscopic origin of this 
doped insulating state, we show in Figs.|3](a) and (d) the 
orbital-resolved DOS for a Na doping of a;=l/8. Close to 
the Fermi level, both dx2-y2 and dxz orbitals contribute 
to the DOS. This can be understood by performing a 
rough energy estimation in the atomic limit when one 
additional electron is added to the Ti t2g bands. While 
the change of total energy is C/ -I- 2 Jh if the additional 



electron occupies the dx2_y2 orbital, it is J7 — Jh -l~ A if the 
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electron occupies the dxz orbital where A is the orbital 
excitation from dx2-y2 to dxz-, and Jh is the Hund's rule 
coupling. Since A ~ 0.1 eV and Jh '~ 1 eV, adding the 
additional electron to the d^z orbital is always preferable. 
This analysis indicates the importance of the Hund's rule 
coupling in the doping case. 

The presented DFT results point to a multi-band sce- 
nario. With the help of A^-th order muffin tin orbital 
(NMTO) downfolding we found that the hopping in- 
tegrals within dxz orbitals along a and within one TiO 
bilayer are comparable to the hopping integrals within 
dx2^y2 orbitals along 6, indicating that the interaction 
network is two-dimensional under doping. Since the ex- 
change interactions among local spins along a and within 
one TiO bilayer are weakly ferromagnetic [13| , the doped 
electron would in principle hop freely in the system 28 1 
and would result in a metallic state. This is though in 
contradiction with the insulating state obtained both in 
our DFT calculations and in experimental observations. 

Therefore, to further elucidate the possible mechanism 
for the insulating state upon doping, we compare the 
crystal field splitting of the undoped case with that of the 
1/8-Na-doped case as shown in Figs.|4](a) and (b). In the 
undoped case all 8 Ti ions in the supercell are equivalent 
and the orbital excitations are the same. Here, the Ti 
ions are labelled from Til to Ti8 as shown in Fig.[T](b). 
In the 1/8-Na-doped case, all Ti ions become inequiva- 
lent and the crystal field splittings are different from site 
to site. Most importantly, the splittings of t2g orbitals on 
the Ti2 (closest Ti to Na) become strikingly small, which 
makes it possible to trap the additional electron into dxz 
orbitals on this Ti site. This finding is consistent with 
the partial DOS shown in Figs. [3] (a) and (d) where the 
peaks near the Fermi level are mainly contributions from 
the dx2-y2 and dxz orbitals of Ti2. The peak at lower en- 
ergies is from the rest of Ti ions (all denoted by Ti2) and 
is of purely (i^2_j,2 character. Further investigations of 
the crystal field splittings and partial DOS for the other 
Na doping concentrations reveal that the same mecha- 
nism can be applied to the observed insulating states, 
i.e., each Na ion strongly modifies the lattice structure 
locally, reduces the crystal field splitting of the t2g or- 
bitals on the closest Ti ion and leads to localization of 
the doped electron, which reduces the Ti^+ {3d^ ) to Ti^+ 
(3(i^) and prevents conduction. Furthermore, long-range 
Coulomb repulsion becomes effective due to the appear- 
ance of mixed- valence state of Ti'^+ and Ti^+ and further 
stabilizes the symmetry-breaking insulator. Within the 
point-charge approximation |16| , the combination of crys- 
tal field splitting and the long-range Coulomb interaction 
can account for the gap of leV. 

We have also considered the effect of substituting O by 
F (Fig. El (b), (e)), CI by S (Fig. [3] (c), (f)), and Ti by 
V or Sc (not shown) in the undoped TiOCl at concen- 
trations of 1/8, leading to TiOy/gFi/gCl, TiOClr/gSi/g, 
TiyygVi/gOCl and TiyygScj^ygOCl. The atomic positions 
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FIG. 4: (Color online) crystal field splitting energies ob- 
tained from GGA calculations. The horizontal axis denotes 
8 Ti ions in the supercell labelled by numbers from 1 to 8 
as shown in Fig. □ (b). (a) TiOCl, (b) Nai/gTiOCl, (c) 
TiOr/gFi/gCl, (d) TiOCl7/8Si/8, (e) Tir/sVi/gOCl, and (f) 
Ti7/gSci/gOCl. 



where the substitutions are made (without loss of gener- 
ality) are specified in Fig. [T] (a). In all the doped cases 
gaps still open at the Fermi level. The DOS of the F- 
doped case (Fig. [3] (b), (e)) is similar to that of the Na- 
doped case with a double peak structure below the Fermi 
level, since both are electron doping processes. The simi- 
larities can be also observed in the crystal field splittings 
(Fig. 131(c)). However, in the F-doped case the Ti2 dx2-y2 
and dxz orbitals are almost degenerate. 

Substituting CI by S implies taking out an electron 
from the dx2_y2 orbital (hole doping), and one would ex- 
pect the system to be metallic. However, from Fig.|4](d), 
we find that the on-site orbital energies of Ti3 (Fig.|4](d)) 
which is closest to the substituted S ion are significantly 
raised due to the large distortion of the lattice structure, 
again leading to a localized state for the doped hole. In 
this case, no spectral weight from the Ti3 3d orbitals is 
detected close to the Fermi level (Fig. [3] (c)). Interest- 
ingly, the Til dx2-y2 orbital energy is almost degenerate 
with the Ti7 dxz orbital energy, indicating the impor- 
tance of on-site Coulomb interaction which avoids double 
occupation on each Ti ion. This leads to a DOS below the 
Fermi level of only dx2-y2 character (Fig. [3] (f)). On the 
other hand, doping with S might result in the formation 
of sir ions, which are isoelectronic to two Cl^ ions and 
do not contribute to doping. The system will then remain 
insulating. Finally, substituting Ti by V (electron dop- 
ing) and Sc (hole doping) we observe a similar effect as 
seen in the previous electron- and hole-doped cases. The 
onsite orbital energies are much lower (higher) on the V 
(Sc) ion than on the Ti ion (Figs.|l](e) and (f)), which can 
again account for the localization of the doped electron 
(hole) and the absence of a metallic state. However, the 
oxidations are aU -1-3 (Ti^+, V^+, Sc^+), indicating that 
the role of long-range Coulomb interaction is negligible. 

From the above analysis, we conclude that the appro- 
priate microscopic model that accounts for the persis- 
tence of the insulating state upon doping is given by a 
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multiband ionic Hubbard model: 



H = ^ tij^rT'c\^„CjT'a + [U + 2 Jh) ^ n„^nir^ 
i.j,a.TT' i,T 



i,r>r',(7 



i,r>r',cr 



(1) 



The first term describes the hopping between Ti sites i,j 
within the three t2g orbitals (r, t'), the second to fourth 
terms are the intra- and inter-orbital Coulomb interac- 
tions and the Hund's rule coupling, respectively, the fifth 
term the crucial on-site orbital energies Ai,r, and the 
last term the effective long-range intersite Coulomb in- 
teraction Vij. Upon Na doping, A^^,- differs from site to 
site which forms a local trapping potential and leads to 
a coexistence of Ti^+ and Ti^+ stabilized by Vij . Then, 
a doping-induced phase transition from undoped Mott 
insulator to doped insulator with charge disproportiona- 
tion occurs. While a similar transition has been exten- 
sively investigated in the one- band case [l3| , inclusion of 
orbital degrees of freedom may result in an even richer 
phase diagram [2^ [s^ . 

Although our DFT calculations are only performed for 
several commensurate Na doping cases, we argue that 
the scenario for this insulating state under doping can be 
applied to all doping concentrations: the Na ion enters an 
individual cage consisting of 5 CI and 1 O ions and creates 
a trapping potential reducing the crystal field splitting of 
the t2g orbitals on the closest Ti ion and resulting in a 
localized state for the doped electron with the help of 
long-range Coulomb interaction. 

Finally, with the understanding of the mechanism for 
the insulating state, we propose two ways to metallize 
TiOCl by doping which would favor a possible supercon- 
ducting state in a doped Mott insulator by suppressing 
other symmetry- breaking states. The first idea is to avoid 
the strong modification of crystal field splitting and the 
formation of a trapping potential due to the lattice de- 
formation induced by the cation or anion. Therefore, a 
possible direction is to intercalate Na together with or- 
ganic ligands to prevent the doped Na from entering the 
cage, in the spirit of intercalating organosolvated Li into 
/3-HfNCl . A second practicable way is to apply exter- 
nal pressure [ij, |15| after carrying out electron doping. 
Due to the appearance of the additional peak and the 
reduction of the gap under doping, the doped system be- 
comes much easier to metallize under pressure than the 
undoped one. 

In conclusion, we performed DFT calculations on 
doped TiOCl and unveiled the microscopic origin of the 
observed insulating state upon doping. The superlat- 
tice structures at different Na doping concentrations and 
the effects of different substitutions are predicted. We 
find that the opportunity to metallize TiOCl or even 



to make it superconducting still remains. Moverover, 
the insulator-to- insulator transition itself has been exten- 
sively discussed by tuning the correlation strength \it\ . 
Here, we open a new way to induce the insulator-insulator 
transition by doping. Finally, our proposed model may 
contain novel phases which calls for further investigations 
by many-body techniques [2^, [s^ . 
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